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ABSTRACT 

We present MMT/Megacam imaging in Sloan g and r of the extremely low luminosity Bootes 
II Milky Way companion. We use a bootstrap approach to perform robust measurements of, and 
uncertainties on, Bootes II's distance, luminosity, size, and morphology. Comparisons with theoretical 
isochrones and empirical globular cluster fiducials show that Bootes IPs stellar population is old and 
metal-poor ([Fe/H] < -2). Assuming a stellar population like that of M92, Bootes II is at a distance of 
42 ± 2 kpc, closer than the initial published estimate of 60 ± 10 kpc. This distance revision, combined 
with a more robust measurement of Bootes II's structure with a Plummer model (exponential model) 
results in a more compact inferred physical half-light size of ~ 36(33) ± 9(10) pc and lower inferred 
luminosity of My ~ — 2.4(— 2.2) ± 0.7(0.7) mag. The revised size and luminosity we calculate move 
Bootes II into a region of size-luminosity space not previously known to be occupied by old stellar 
populations, but also occupied by the recently discovered Milky Way satellites Willman 1 and SEGUE 
1. We show that the apparently distorted morphology of Bootes II is not statistically significant given 
the present data. We use a tidal argument to support a scenario where Bootes II is a dwarf galaxy 
(dark matter dominated) rather than a globular cluster (not dark matter dominated), although the 
uncertainty on the M/L we infer for Bootes II is substantial. Moreover, we can not rule out that 
Bootes II is a star cluster on the verge of disruption, such as Palomar 5. 
Subject headings: galaxies: dwarf — Local Group 



1. INTRODUCTION 

Over the last 5 years, the Sloan Digital Sky Survey 
has been extensively searched for extremely low surface 
brightness dwarf spheroidal galaxies. These searches use 
the catalog of stellar sources to identify spatial overden- 
sities of the old, metal-poor stars c haracteristic of these 
dwarfs. To date, thes e searches fe.g lWillman et aDl2002t 
iKoposov et al.ll2007al ) have resulted in the discoveries of 
fourteen new Milky Way satellites. 

Most of these new objects have total luminosities 
less than the median luminosity of the Milky Way's 
globular clusters (My ~ —7), but have sizes charac- 
teristic of known dwarf spheroidals (r^aif ^ 100 pc), 
complicating their classifications as either star clus- 
ters or dwarf galaxies. In this paper, we assert that 
the physical distinction between a globular cluster and 
a dwarf galaxy is that a dwarf galaxy is, or was 
at some point, the primary baryonic component of a 
dark matter halo whereas a globular cluster was not. 
Nine of the 14 new satellites were originally classi- 
fied as dwarf spheroidals because they had scale sizes 
> 100 pc (Bootes, Canes Venatici, Canes Venatici II, 
Coma Berenices, H ercules, Leo IV, Leo T , Ursa Major 
and U r sa Major II: iBelokurov et al.l 120061: IZucker et al 
2006ri ISakamotcT fc Hascgawa 2006; Belo kurov et al 



appear to require dark matter to explain their kine- 
matics, securing their classification as dwarf galaxies 
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iMartin et all 120071 : IStrigari et all 
2007a; Wal ker et all I2007D. Koposov 1 and 2 
(|Koposov et al.H2007bf ) have scale sizes of only 3 pc, and 
have thus been classified as new globular clusters. 

The classification of the other three new satellites, 
Segue 1, Willman 1, and Bootes II (M v ~ -2.5; 



Belok urov et~alll2007t IWillman et aT1l2005at IWalsh et all 



2007tlWillman et al.l2005btlZucker et alJl2006aD . Follow- 
up spectroscopic studies demonstrated that they indeed 
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20071) . has been less straightforward. They are old stel- 
lar populations with sizes intermediate between known 
Milky Way globular clusters and dwarf spheroidals, but 
have fewer stars than nearly any known galaxy or glob- 
ular cluster. Although initial estimates based purely on 
SDSS data placed Bootes II close to Coma Berenices in 
size-luminosity space at (log(r/t/pc), My) = (1.85,-3.1), 
in this paper we present estimates based on a more ro- 
bust algorithm and on deeper, MMT/MegaCam imaging 
in g and r that shift Bootes II's size and luminosity closer 
those of Willman 1 and Segue 1 . 

Despite its tiny luminosity, spectroscopic [Fe/H] es- 
timates and kinematic studies of Willman 1 have pro- 
vided support that it may require dark matter to ex- 
plain its properties , thus classifying it as a dwarf galaxy 
(|Martin et aU2007tlStrIgari et al.ll2008f) . Whether or not 
these three Milky Way companions with only ~ lOOOi© 
are galaxies or globular clusters is of fundamental import 
to both our understanding of galaxy formation at the 
smallest scales and to our understanding of the size and 
mass scale of dark matter clustering. Although the ex- 
tent to which tides have effected Willman l's present-day 
luminosity is uncertain, this object raises the questions: 
Are we for the first time seeing the low luminosity and 
mass limit of galaxy formation? If so, what properties 
can wc infer for the dark matter halos that host such 



2 



Walsh et al. 



galaxies? Segue I and Bootes II presently lack their own 
published spectroscopic studies with which to evaluate 
these scenarios. However, rigorous derivations of the de- 
tection l imits of the most rece nt SDSS searches for such 
objects (|Koposov et al.ll2007at Walsh, Willman & Jerjen 
in prep) show that many similar objects could remain 
yet undiscovered around the Galaxy, underscoring the 
importance of understanding their physical properties. 

With the study presented in this paper, we aim to pro- 
vide the first robust measurements of the basic properties 
(distance, luminosity, structure) of Bootes II and its stel- 
lar population, and to evaluate the present evidence for 
its classification. In §2, we describe MMT/McgaCam ob- 
servations of Bootes II, data reduction, and artificial star 
tests. In §3, we use these data to derive revised estimates 
of Bootes IPs properties, to verify that its stellar popu- 
lation is old and metal-poor, and to investigate whether 
it has a distorted morphology. We discuss these results 
and evaluate evidence for a dwarf galaxy versus globular 
cluster classification of Bootes II in §4. 

2. DATA 

We observed Bootes II on June 05 2007 with Mega- 
Cam (jMcLeod et al.ll2000t ) on the MMT. These data were 
obtained as part of a larger survey program to image 
with MMT/Megacam ultra- faint Milky Way satellites. 
MMT/MegaCam has 36 chips with 2048x4608 pixels of 
0.08"/pixel, for a total field-of-view (FOV) of 24'. We 
obtained 5 180s dithered exposures in Sloan g, and 5 
240s dithered exposures in Sloan r in grey conditions 
with 1.0 - 1.2" image quality in the g and 0.9 - 1.0" 
image quality in the r images. We reduced the data 
based on the method described in Matt Ashby's Mega- 
cam Reduction Guide 1 . Our reduction relied in part on 
software written specifically for MMT/MegaCam data 
reduction by Brian McLeod. We used the Sloan Digital 
Sky Survey Data Release 6 (SPS S DR6: lYork etliI1l2000l : 
lAdelman-McCarthv et a l. 2008) stellar catalog to derive 
precise astrometric solutions for each science exposure. 
We also used the SDSS catalog to derive an illumination 
correction in g and r to divide out the variation in zero- 
point across MegaCam's FOV. We use local copies of the 
SDSS dataset, maintained at the Harvard-Smithsonian 
Center for Astrophysics. 

We did a weighted co-addition of the reduced images 
using S WARP 2 and th en used the DAOPHOTII/Allstar 
package (|Stetsonlll994h to do point source photometry on 
the resulting images. We visually verified the integrity 
of the shape and full-width half-max of the PSFs in the 
stacked images across the 24' FOV. Photo metry was car - 
ried out using a method similar to that of lHarrisl (|2007h . 
with the exception that we used the command-line ver- 
sions of DAOPHOT and Allstar rather than the IRAF 
versions. 

To derive the photometric calibration for our data, we 
first matched the SDSS stellar catalog to the Allstar cat- 
alog for these new observations. We used the 91 SDSS 
stars within our field-of-view (FOV) with 18 < r < 21 
and 0.1 < g — r < 0.8 to perform the photometric cali- 
bration. We limited the calibration to stars fainter than 
r = 18 mag to avoid the saturation limit of the Mega- 

1 http:/ /www. cfa.harvard.edu/~mashby/mcgacam/mcgacam_frE 

2 http:/ /terapix. iap.fr/soft/swarp 



Cam data. We limited the calibration to stars with col- 
ors between 0.1 < g — r < 0.8 because the Bootes II 
member stars resolved in this study (with the exception 
of a few possible blue horizontal branch stars) all have 
0.1 < g — r < 0.6. There were insufficient SDSS stars in 
our FOV bluer than 0.2 mag to determine whether our 
derived calibration is appropriate for very blue stars. 

We then did a linear least-squares fit for the zero- 
points and color-terms, including uncertainties in color 
and magnitudes on each star and throwing out 3 sigma 
outliers. 

9 = 9instr + 7.27(±0.029) + 0.091(±0.068) x (g - r) (1) 

r = r instr + 7.33(±0.025) + 0.074(±0.054) x(g-r) (2) 

Uncertainties were derived from a 1000 iteration boot- 
strap of the data. In addition, there is uncertainty 
in the SDSS zero-points th emselves of about 0.01 mag 
(jPadmanabhan et al.ll200cl . 

Throughout this paper, we adopt SDSS photometry, 
rather than MegaCam photometry, for stars brighter 
than r = 18.0 mag. All magnitudes in this paper 
have been ext i nction corrected with the values from the 
ISchleeel et~aH (|1998l ) dust maps provided in the SDSS 
catalog; the median E(g — r) along the line-of-sight to 
Bootes II is 0.02 mag. 

We use artificial star tests to measure the photometric 
errors and completeness as a function of position in the 
g — r color-magnitude diagram (CMD). Artificial stars 
are constructed from the g and r point spread functions 
(PSFs) measured during the data-reduction process, and 
are injected into the co-added g and r images using a uni- 
form grid with spacing in X and Y equal to ten times the 
full-width half-max (FWHM) of the PSF (so that artifi- 
cial stars overlap only beyond their 10cr radii). This fixed 
geometry imposes a limit on the number of artificial stars 
that can be added to the image of about 18,500. To build 
up our number statistics, we inject artificial stars into 
twenty copies of the g and r images, randomly offsetting 
the grid's zero-point position in X and Y for each itera- 
tion. This results in a total sample of 370,000 artificial 
stars. The r photometry of the artificial stars is drawn 
randomly from ~ 18 to 28 mag, with an exponentially in- 
creasing probability toward fainter magnitudes. To prop- 
erly characterize the tail of the completeness function and 
the impact of blends on the photometric errors of faint 
objects, we simulate stars up to three magnitudes fainter 
than the nominal faint limit. The g — r color is then 
drawn randomly over the range -0.5 to 1.5 mag to deter- 
mine the g magnitude. We photometer the artificial-star 
images with the same photometry pipeline as we used on 
the science frames. 

If an input artificial star is not present in both the g 
and r Allstar files, then it is flagged as a non-detection 
for calculation of the completeness rate. We applied a 
strict cut of DAOPHOT sharpness parameter of — 1 < 
sharp < 1 for a star to be included our analysis, both 
for the actual Boo II data and for the artificial star tests. 
Although this strict sharpness cut yields completeness 
limits that are brighter than if we use no cut at all, we 
found it necessary to eliminate many galaxy interlopers 
Hind provide an improved measurement of Bootes IPs stel- 
lar population while not sacrificing much precision in our 
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quantitative results. At a g — r color of 0.25 mag, the 
50% and 90% completeness limits of these data are r = 
23.5 and 22.9 mag. 

Figure [T] shows the CMDs of stars within 4.4 arcmin- 
utes of the center of Bootes II from the SDSS DR6 data 
and from the MMT/Megacam data. The 50% complete- 
ness, as a function of color, is overplotted on the Mega- 
cam CMD. 

3. BOOTES II PROPERTIES 

3.1. Bootstrap Analysis 

With only ^100 object stars resolved in this study, 
small number statistics will constitute a substantial, if 
not dominant, source of uncertainty in the derived quan- 
tities. For the ambiguous ultra-faint satellites such as 
Bootes II, a rigorous examination of the uncertainties 
is essential to measure any of their properties. We use 
a 10,000 iteration bootstrap analysis to determine both 
the most likely values of Bootes II's properties, and the 
uncertainties associated with each measurement as de- 
tailed in the following sections. For each iteration, the 
data are randomly re-sampled with replacement and then 
analyzed to derive as described in the remainder of §3: 
central RA and Dec, distance modulus, Plummer and ex- 
ponential half-light radii and total absolute magnitudes, 
King core and tidal radius, position angle, cllipticity, 
and asymmetry Aside from King tidal radius, the boot- 
strapped distributions of the derived parameters are well 
described by a Gaussian. All quoted values and uncer- 
tainties are thus the peaks and standard deviations of 
the bootstrap distributions. For the King tidal radius, 
we quote the half-width, half-max as the uncertainty be- 
cause there is a long, poorly populated tail of values that 
extends to high tidal radius. We report the values of 
bootstrapped sample fits in Table 1. 

3.2. Central Position 

Bootes II contains no detectable unresolved luminous 
component, so we determine its center by locating the 
barycenter of likely member stars within a 4.5 arcmin ra- 
dius around the center published in Walsh et al. (2007). 
We define likely member stars as those that have col- 
ors and magnitudes consistent with M92 at an initial 
estimated distance modulus. We repeat this with the 
returned RA and Dec until the difference between input 
and output values converges. Our derived values are pre- 
sented in Table 1 along with their associated bootstrap 
uncertainty. The uncertainties on the center of Bootes II 
are substantial: 7.2 seconds of RA and 23 arcseconds in 
Dec. Because all parameters are derived for each of the 
bootstrapped samples, these uncertainties on the center 
are automatically propagated through to the uncertain- 
ties in Bootes IPs structural parameters. 

3.3. Distance 

To investigate the distance to Bootes IPs stellar pop- 
ulation, we first compare its CMD to empirical glob- 
ular cluster fiducials (M92, M3, M13, and M71) with 
-2.4 < [Fe/H] < -0.7. We use m - M = 14.60, 15.14 , 
14.42 and 13.71 for t he four clusters dPaust et iD|2007|; 
Kraft fc Ivans! 120031 iCho et all [20051 iGrundahl et all 
20021 ). We choose to rely on fiducials, rather than the- 
oretical isochrones, because these well studied globular 



clusters have photometry in the exact photometric sys- 
tem we have calibrat ed our data to. Th e fiducials we use 
are based on those of IClem et al.l (|2008h in Sloan g' — r' . 
They were converted int o g — r using the transforma- 
tion of lRider et al.l (|2004l ) and checked by comparing the 
transformed fiducials directly to the SDSS imaging of of 
the clusters in the SDSS DR6 (J. Strader, private com- 
munication) . The robustness of our comparison depends 
on Bootes II having an old stellar population, like those 
in these four comparison clusters. We address and con- 
firm this with theoretical isochrones in §3.4. 

For each fiducial, we find the distance modulus that 
provides the best fit to the stars in the CMD shown in 
the right panel of Figure [1] This CMD includes all stars 
within 4.5 arcmin of Bootes II. To determine this distance 
modulus, we step each fiducial through 0.05 magnitude 
intervals in (m — M) from 17.5 to 19.5 mag and find the 
number of stars brighter than r = 23.5 that, considering 
color uncertainties, have colors within 0.05 magnitudes 
of the fiducial. To eliminate the contribution of stars 
belonging to the thick disk and halo, we then do the same 
for stars exterior to 9.0 arcmin and subtract this value 
normalized to an area of 7r(4.5') 2 . We take the best fit 
for each fiducial as the distance modulus that maximizes 
this number of stars. The best fit distance moduli for 
the M92, M3, M13 and M71 fiducials are 18.1, 18.1, 18.1 
and 18.85 with 96, 77, 87 and 43 stars respectively. 

The left panel of Figure [U shows the M92, M13 and 
M71 fiducials overplotted on the MMT CMD at their in- 
dividual best fit distances. For simplicity we include M13 
and not M3, because they have similar [Fe/H] and M13 
provides a better match both quantitatively and quali- 
tatively. In the central panel of Figure [2] we overplot the 
empirical M92, M13 and M71 fiducials, all projected to 
the M92 best fit distance modulus of m-M = 18.1 (42 
kpc). These confirm that it is reasonable to infer that 
the stellar population of Bootes II is like that of M92. We 
therefore use the M92 fiducial in our bootstrap analysis 
to derive the best-fit distance, which yields a distance 
modulus of m — M = 18.1 ± 0.06, only including the for- 
mal bootstrap uncertainty, stemming from small number 
statistics. We add in quadratur e the distance mo dulus 
uncertainty of M92 (0.09 mag, iPaust et all l2007f ) and 
the r zero-point uncertainty (0.025 mag, §2) to derive 
m-M = 18.1 ± 0.06, or d = 42 ± 2 kpc. If Bootes 
II has a stellar population different from M92, then the 
uncertainty in the distance is larger. 

3.4. Stellar Population 

M92 has a very low [Fe/H] of -2.4 and is a-element 
enha nced relative to solar ([Ca/Fe] = 0.3. ISneden et all 
I2000T ). The match between M92 and Bootes IPs stel- 
lar populations thus supports [Fe/H]s 00 // < —2, even if 
Bootes II is a-depleted relative to M92 (typical of the 
contrast between dSph and globular cluster populations, 
iPritzl et al.1l2005D . Figure [I shows that M13's fiducial 
sequence is only 0.04 mag redder in (g — r) than that of 
M92 below the main-sequence turnoff. The uncertain- 
ties in the g and r zero-points of the MegaCam data 
in Figure O result in a (g — r) calibration that is un- 
certain at 0.038 mag. These Bootes II data are thus 
also consistent with having a more moderate abundance 
{[Fe/H] ~ —1.6). However, using independent V and I 
observations obtained on VLT/FORS2, Jerjen et al. (in 
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Fig. 1. — Color-Magnitude Diagrams of stars within 4.5 arcmin (~ 1.5 X Plummer half-light radii) of the center of Bootes II based on 
MMT/Megacam observations. The left CMD shows data from the SDSS DR6 and the right panel shows the MMT/Megacam data. Error 
bars showing the color and magnitude uncertainties as a function of r are overplotted. 
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Fig. 2. — Color-Magnitude Diagrams of stars within 4.5 arcmin of the center of Bootes II Left: Globular cluster fiducials overplotted 
at their own best fit distance modulus. Center. GC fiducials overplotted at M92 best fit distance modulus of m — M = 18.1 Right: M92 
fiducial at m — M = 18.1 with probable Bootes II member stars highlighted. 



preparation) also find that Bootes II is best described by 
the most metal-poor fiducials. 

To check these empirical results and to investigate a 
range of possible ages for Bootes II, we repeat the dis- 
tance modulus fitti ng described i n §3.3 using the theo- 
retical isochrones of iDotter et al.l (|2008f ) in SDSS colors. 
We use 24 isochrones corresponding to the combinations 
with [Fe/H]= -2.3, -1.5 and -0.7, ages of 5, 7, 10 and 
13 Gyr and alpha-element abundances of [a/Fe] = 0.0 
and 0.2. The best fitting of all 24 isochrones is that with 
[Fe/H]= -2.3, 13 Gyr and [a/Fe] = 0.2, with 98 stars 
having colors lying within 0.05 mags of the isochrone. For 
the same abundance values, the number of stars drops 
to 86, 78 and 72 for the 10, 7 and 5 Gyr populations re- 
spectively. This quantitative comparison with the Dotter 
isochrones highlights the fact that the small color differ- 
ence between the MSTO and RGB stars in Bootes II 



would not be consistent with a stellar population much 
younger than 13 Gyr. 

3.5. Structural Parameters 

The surface density profiles of globular clusters and 
dwarf spheroidal galaxies (dSphs) are co mmonly param- 
eterized by King (|Kinglll966l ). Plummer (|Plummerl 1911) 
and exponential profiles. To facilitate comparison with 
other observational studies, we fit all three profiles to the 
stellar distribution of Bootes II: 
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where rp and a are the scalelengths for the Plummcr 
and exponential profiles and r c and r t are the King core 
and tidal radii, respectively. For the Plummer profile, 
rp equals the half-light radius rn, while for the expo- 
nential profile r/j ~ 1.668a. The circled stars in the 
right panel of Figure [5] show the color-magnitude crite- 
ria we use to select probable Bootes II member stars for 
calculating its center and for investigating its structure. 
Figure [3] shows the spatial distribution of stars that pass 
these cuts, and the location of their derived center. Fig- 
ure[4]shows the surface density profile of Bootes II around 
this center, where the error bars were derived assum- 
ing Poisson statistics. Using a non-linear least squares 
method, we fit Plummer and exponential models plus a 
constant field contamination to this surface density pro- 
file. The surface density profile fits are only constrained 
to be physically possible systems (i.e. field contributions 
must be positive). In the case of the King profile, the 
tidal radius is present in a constant term, hence there 
is a degeneracy between the field value and the tidal ra- 
dius. To circumvent this we first fix the King field value 
using the mean of the fitted Plummer and exponential 
field values. 

All three fits yield consistent characteristic radii; As- 
suming a distance modulus of m — M — 18.1, the Plum- 
mer and exponential profiles yield physical half-light radii 
oi r h ,pi urnmer ~ 36 ±9 pc and ^exponential - 33 ±10 pc. 
The King model fit yields a core radius of r c ~ 25 ± 9 pc 
and a tidal radius of r t ~ 155±35 pc. Although this tidal 
radius lies just outside the extent of the radial profile we 
can measure, we find that the inner radial bins constrain 
the core radius and central density while fixing the field 
value leaves only the tidal radius as a free parameter. 
In the event that the outer radial bins are contaminated 
with Boo II stars, and therefore higher than the true field 
value, the tidal radius could be much larger. As an ex- 
ample, if the true field value is an overestimate of 25% by 
the Plummer and exponential models, then the best fit 
tidal radius is ~ 215 pc while the core radius remains rel- 
atively constant at ~ 22pc. The fitted King tidal radius 
of r t ~ 155 ± 35 is a lower limit. 

Figure [5] shows in greyscale the distribution of King 
core and tidal radii derived from the bootstrap of Bootes 
II stars. Overplotted are the measured core and tidal 
radii of Galactic globular clusters. We have overplot- 
ted lines of constant concentration (rudai/r C ore) and the 
concentrations of known globular clusters calculated wit h 
the tidal and core radii in the catalog of iHarria fl996). 
Only a handful of known globulars have > 10 pc and 
the only one larger than Th ~ 20 pc is known to be 
tidally disrupting (Pal 5). However, this figure shows 
several known globular clusters with a King concentra- 
tion as low as Bootes II. The King concentrations of the 
Milky Way's dwarf spheroidals range from 0.48 - 1.12 
(|Mateolll998f L similar to the range of GCs in this figure. 
We do not overplot the King parameters of the classical 
Milky Way dwarf galaxies, because their relaxation times 
are ~ a Hubble time. Their King profile fits thus contain 
less physical meaning than those of objects with shorter 
relaxation times, such as globular clusters and Bootes II 
(see §4 for discussion). 



-0.1 - 



-0.2 t_ 




0.1 0.0 -0.1 -0.2 

Aa (deg) 

Fig. 3. — Positions of stars passing selection cuts. The calculated 
center of Bootes II is highlighted by the crosshair, which spans the 
RA and Dec uncertainties. The large circle shows the maximum 
radius of the surface density profile fit. 
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Fig. 4. — Fitted surface density profile of Bootes II consisting 
of a Plummer (blue), Exponential (red) or King (black) profile 
combined with a constant field contribution. 

3.6. Luminosity 

To estimate the total luminosity of Bootes II we inte- 
grate the model components of the surface density pro- 
file, corrected for incompleteness, to derive the number of 
Bootes II stars within the CMD selection cuts. Using this 
to normalize the theoretical luminosity function gives es- 
timates of M r = —2.6 (My = —2.4) for the Plummcr 
profile and M r = —2.4 (My = —2.2) for t he exponential 
profil e (using V — r = 0.16, adapted from lGirardi et aTl 
120041 for a 13 Gyr, [Fe/H] = -2.27 stellar population) 
. We obtain these magnitudes after correcting for the 
missing flux from stars fainter than r = 23.5 by inte- 
grating the theoret ical luminosity function taken from 
iGirardi et all IpOOl for a [Fe/H]= -2.27, 13 Gyr popu- 
lation. This method counts the number of stars without 
regard to their individual magnitudes, which in systems 
of such low luminosity could strongly be affected by the 
addition or subtraction of a single RGB star. Such a 
star's individual magnitude could be much brighter than 
My = —2, rendering a more traditional summing of 
fluxes method unreliable. As the total luminosity of a 
stellar system becomes comparable to the luminosity of 
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Fig. 5. — 2D Histogram of bootstrap results for Bootes II's 
King core radius and tidal radius with Galactic globular clusters 
(GGCs). GGCs having r c > lOpc or rt > lOOpc are shown with 
larger points. Dashed lines show contours of constant concentra- 
tion, defined as c = log(rt/r c ). The open circles show the 8 glob- 
ulars with 35 < d < 100 kpc for which we calculate M/L in §4.3. 

individual stars, the summed luminosity becomes dom- 
inated by the brightest stars which may or may not be 
present simply from small number statistics. Hence, a 
group of stellar systems with an identical number of stars 
in each will have a spread in luminosity independent 
of any measurement uncertainty (see also iMartin et al.l 
(2008)). Using simulated dwarf galaxies from Walsh, 
Willman & Jerjen [in prep) we find that for an object 
such as Bootes II as observed by MMT/Megacam, this 
spread has a standard deviation of 0.6 magnitudes. We 
therefore combine this effect with the bootstrap uncer- 
tainty by summing in quadrature to derive the values 
presented in Table 1. 

3.7. Morphology 

We look for evidence of tidal disturbance based on the 
morphology of Bootes II's isodensity contours. Binning 
the positions of stars in Figure [3] into 0.01° x 0.01° bins 
and spatially smoothing with a Gaussian of 1.5 arcmin 
FWHM scale length reveals an apparent distorted mor- 
phology to Bootes II, including substructure at the 3-5 
a level and an elongation directed along the gradient of 
the Galactic potential, as shown in Figure [HI Although 
tidal tails trace an object's orbit, tidal debris is stripped 
from an object along the gradient of the gravitational po- 
tential such that tidal stars near an object are expected 
to lie along this gradient. The elongation of Bootes II 
could thus be a feature resulting from tidal interaction. 

However, due to the meager number of stars in Boo II, 
any observed irregular morphology may be an effect of 
small number statistics. To evaluate the significance of 
the morphology shown in Figure [6j we generate contour 
plots of bootstrap resamples of Bootes II stars. Figure 
[7] shows nine such randomly selected isodensity contours 
defined the same way as in Figure [51 This figure shows 
that the irregular morphology and apparent distortion 
along the Galactic potential are not persistent features. 
While this does not rule out tidal disturbance to Bootes 
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Fig. 6. — Smoothed contour plot of Bootes II showing an ap- 
parent distortion along the direction towards the Galactic Center. 
The contours show the 3,4,5 and 6 sigma levels above the mean. 

II, the varying morphologies of these resampled objects 
demonstrates that the shape of Bootes II does not neces- 
sarily reflect a true irregularity in its underlying spatial 
distribution. 

In order to quantify any asymmetry in the morphology 
of Bootes II, we first derive the position angle 9 and ellip- 
ticity e. W e calculate these from the standard SExtractor 
definitions (|Bertin fc Arnoutil 19961 ), using the smoothed 
images assuming that all pixels greater than 3tr above the 
mean are part of Bootes II. We then count the number of 
stars on either side of the major and minor axes, within 
1.5r/j. If the positions of these stars are drawn from an 
axisymmetric distribution, then the numbers on cither 
side of an axis should be within y/2(N) of each other. 
We define an asymmetry parameter A: 



A = 



N, - No 



(G) 



where Ni and N2 are the counts on either side of the 
axis and (N) is their average. Hence, for each bootstrap 
iteration we have two values of A, one for the major 
and one for the minor axis. Doing this for simulated 
Boo Il-like objects drawn from a pure Plummer den- 
sity profile expectedly yields a distribution of A with a 
mean of 0.0 and a standard deviation of cr sirns = 1.01. 
The bootstrap yields two distributions with means of 
Aminor = —0.23 and A ma j or — —0.18 and standard devi- 
ations of o m inor = 1-03 and o- ma jor = 1.01. The asymme- 
try of Bootes II is therefore not statistically significant. 

These results imply that any apparent asymmetry in 
the distribution of Bootes II stars is well within that ex- 
pected from a symmetric system, and is probably due to 
small number statistics. Deeper and wider-field imaging 
may provide the signal necessary to definitively measure 
whether Bootes II has extended, tidal structure. 

4. IS BOOTES II DARK MATTER DOMINATED? 

As discussed in the Introduction, one motivation for 
studying an individual ultra-low luminosity object such 
as Bootes II in great detail is to determine whether 
it is a galaxy or a star cluster. Using our definition 
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Fig. 7. — Smoothed contour plot of 9 randomly selected resamples of Bootes II stars, showing varying morphology for each iteration. 
The contours show the 3,4,5,6 and 7 sigma levels above the mean. 



of a dwarf spheroidal galaxy, this boils down to deter- 
mining whether or not the object is dark-matter domi- 
nated. (However, for exam ple, see lMetz fc Kroupall2007t 
Dabringha usen et al.l [200 8 for another interpretation of 
the Milky Way's dwarf companions). Understanding the 
properties of dark matter depends critically on identify- 
ing the smallest ma ss and length scales at which dark 
matt er clusters (e.g. IStrigari et al.ll2007bt iGilmore et all 
l2007Tj . The most direct way we have at present to inves- 
tigate dark matter on the smallest scales is by using the 
least luminous galaxies as tracers of dark matter, and by 
pushing the envelope to find the smallest galaxies possi- 
ble to form. 

In this section, we use a tidal argument to show why 
Bootes II may be dark matter dominated, and thus 
a dwarf galaxy, despite its very low luminosity (L ~ 
500Lq). The lines of evidence presented in this section 
are circumstantial now, but are illustrative of arguments 
that could provide strong constraints on the correct clas- 
sification for Bootes II, if deeper, wider-field imaging, 
kinematics, and/or spectroscopic abundance measure- 
ments are able to demonstrate whether Bootes II is self- 
bound or not. 



4.1. Separation in size-luminosity space from dwarf 
galaxies and globular clusters 

In Figure [51 we compare the size and luminosity 
of Bootes II with that of other nearby, old, stel- 
lar populations. This figure shows Bootes II in the 
size-lumi nosity plane along with Milky Wa y globular 
cluste rs ([Harris! 119961: IKoposov et ail l2007bl. GC s) and 
M31 (|McConnachie & Irwinl 120061 : IZucker et all 2004 
20071) and Milky Way dwarf satellites (iMateol 11991 
Grebel et al.ll2003t iBelokurov et alj|2007t IWillman et al l 
2005bt IZucker et all I2006al: Willman et al.l l2005aF 
Irwin et all 120071 : IZucker et all l2006bl: IBelokurov et all 
2006h . In this figure, including only known old stellar 
populations, Bootes II, Will and Segue 1 occupy a 
somewhat unique place in size-luminosity space. The 
apparent lack of objects with half-light radii larger than 
those of Bootes II, Wl and Segue 1 but smaller than the 
apparent minimum size of confirmed dwarf galaxies (~ 
100 pc) is quite possibly an observational selection effect 
due to the faint central surface brightnesses of objects 
in that region of the size-luminosity plane. 

Conversely, Bootes II, Will, and Segue 1 have half-light 
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Fig. 8. — Size-luminosity plot of known Galactic satellites. The 
previous position of Bootes II is shown in grey along with the region 
bounded by the uncertainties in r h and My- Values f or the new 
Milky Way satellites are taken fromlMartin et al. (2008). Globular 
cluster properties are from Harris (1996). Dotted lines show lines 
of constant surface brightness. 

radii (20 - 40 pc) an order of magnitude larger than the 
half-light size characteristic of similarly low luminosity 
globular clusters (—4 < My < — 1.) Low luminosity 
Milky Way clusters in that siz e gap could have been de- 
tected by the SDSS searches of lKoposov et al.1 (|2007aD or 
Walsh, Willman & Jerjen (in prep). This gap between 
Bootes II, Willman 1 and Segue 1 and the known globular 
clusters is thus real and not a selection effect, suggesting 
that these objects are a distinct population from Milky 
Way GCs. However, if such diffuse old clusters do exist 
outside of the Milky Way, they would not yet have been 
discovered owing to their low surface brightnesses. It 
thus remains possible that old star clusters with proper- 
ties bridging the gap between Bootes II and known GCs 

do exist in abundance in other e nvironment s. 

Based on the compilation of iDias et al.1 (|2002l ) (cat - 
alog obtained at |http: / /www.astro.iag.usp.br~wilton/| ) , 
only 2 of the 1076 Milky Way open clusters in that 
catalog with both distance and diameter measurements 
have apparent radii larger than 20 pc. Young stellar 
associations in environments other than the Milky Way 
have been observed to have characteristic sizes of up to 
100 pc (e.g. OB asso ciations in the LMC observed by 
iGouliermis et al.ll2003f ). Although the majority of these 
clusters appear to be unbound, some fraction of them 
could survive as self-bound entities for 10 Gyr if subject 
only to very weak tidal forces. We will show in §4.2 and 
4.3, that tidal effects limit the survivability of objects at 
Bootes II distances with star cluster-like mass-to-light ra- 
tios and sizes and luminosities similar to those of Bootes 
II, Willman 1, and Segue 1. 

4.2. Minimum mass-to-light ratio to be bound 

We estimate a lower limit to the mass-to-light ratio 
required for Bootes II to be self-bound. We calculate the 
instantaneous tidal radius of Bootes II using: 



rudai a 



/ rriBooii \ 
\3M M w(d)J 



>*BooII, 



(7) 



where rasooii is the mass of Bootes II, MMw{d) is the 
mass of the Milky Way within the distance to Bootes II, 
and oIbooII is the galactocentric distance to Bootes II. 
We calculate the Milky Way mass Mmw (d) assuming an 
isother mal sphere with a circular velocity of V c = 220±40 
km s^ 1 (|Bellazzinill2004l ) . using dsooii = 42 ±2 kpc, and 
using M v = -2.4 ± 0.7 mag for Bootes II . For (M/L) v 
values of 2, 10, 100 and 1000 kpc the corresponding tidal 
radii are thus 42±^, 72±^, 156±^ and 336lg^ pc re- 
spectively. A tidal radius of 72 pc (M/L = 10) is within 
our observed field and only ~ 2 times the half-light radius 
of Bootes II. The visible extent of Bootes II thus exceeds 
that expected for globular cluster-like mass-to-light ra- 
tios. If the Bootes II stars throughout our field-of-view 
are indeed bound to the object, then it requires a signif- 
icant dark matter component. In §3.7, we showed that 
Bootes II lacks statistically significant tidal distortion, so 
there is presently no evidenc e that is is losing stars (how- 
ever see lMunoz et al.l 1200 8). Deeper wide- field imaging 
and/or kinematic data will provide a clearer picture of 
the boundedness of this object. 

4.3. Mass-to-light ratio inferred from the King tidal 

radius 

Using the same physical principle as in the previ- 
ous section, we now calculate the mass-to-light ratio of 
Bootes II assuming that the King tidal radius derived in 
§3.5 is the true tidal radius of Bootes II, ow ing to the ti dal 
field of the Milky Way. The King model (|Kinglll966h is 
physically motivated and expected for a relaxed, single 
mass component, spherical system in equilibrium that 
is tidally limited by the gravitational field of the Milky 
Way. For velocity dispersions greater than 0.1 km/sec, 
the relaxation time of Bootes II is less than a Hubble 
ti me (unlike th e class ical Milky Way dSphs). Based on 
a iDotter et al.l (|2008f ) stellar luminosity function for a 
13 Gyr, [Fe/H] = — 2 system normalized to have My 
— -2.4 mag, Bootes II has 3700 stars more massive than 
0.1M,q„,Ti.. Its relaxation ti me (t re i ax ~ N/8ln(N)*t cross ; 
iBinnev fc Tremaineill987l ) at its Plummer half-light ra- 
dius of 36 pc is thus ~ 1.1 Gyr/<x„, where a is its ID 
velocity dispersion. Although the presence of a signif- 
icant dark component of matter would complicate the 
expected relaxation time of Bootes II, we are testing the 
hypothesis that it does not contain dark matter. 

If we take the King "tidal radius" to be the instan- 
taneous tidal radius of Bootes II, we can infer a mass 
mBooii and therefore a mass-to-light ratio. We use Equa- 
tion 7 as described in §4.2, but set r tidal = r tidaliKlng = 
155 ± 35 pc, and solve for the mass of Bootes II. We infer 
a V-band mass-to- light ratio of Bootes II of 98^4°. If 
this is representative of the true Mj L of Bootes II, its 
mass would be dominated by dark matter, classifying it 
as a galaxy. We note the substantial uncertainty on the 
M/L derived in this way. Even if the input assumptions 
are robust, the dark matter signal is only significant at 
the 1-tr level. As mentioned in §3.5, the King tidal radius 
may be substantially larger if the outer surface density 
bins are contaminated with Boo II stars thereby overes- 
timating the field density. As in §3.5, if the field value 
is overestimated by 25%, then the best fit tidal radius 
would be r t = 215 pc and the inferred M/L would be 
263+ 991 
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For comparison, we use these same assumptions and 
same technique to calculate the mass-to-light ratios of 
known Milky Way globular clusters in the distance range 
of 35 < Rgc < 100 kpc. We chose this distance range 
because i) these halo globulars reside in an environment 
and are on orbits that are th e most comparable to Bootes 
II, and ii) iBellazzinil (|2004) states the Milky Way mass 
profile is consistent with an isothermal sphere with a cir- 
cular velocity of V c = 220 ± 40 kms^ 1 between 35 - 100 
kpc. There are 8 clusters in that range: Eridanus, Pal 
2, NGC 2419, Pyxis, Pal 3, Pal 14, Pal 15 and NGC 
7006. For all of these cluster s we u se the distances and 
luminosities from the lHarrisI (|1996l ) catalog, and assign 
a distance uncertainty of 0.1 mag in distance modulus. 
For Pal 2, NGC 2419, Pal 3, Pal 14 and NGC 7006 we 
use the core radii, concentrations an d uncertainties from 
McLaughlin & van dcr Marcl (2005). For the other three 
GCs. we use values from lHarrisI (|1996t ) and assign uncer- 
tainties of 10 pc to their tidal radii. We find a median 
(mean) (M/L) v of 0.36±^ (0.5±J;||) for these eight 
clusters. The single outlier with a calculated M/L > 1 is 
Pal 14 with M/L = 2A±{° 9 4 . Although the true M/L of 
most (all) of these eight halo globulars is within lcr (2a) 
of the calculated values, we find that this technique sys- 
tematically underestimates of the M/L of these relaxed 
systems by up to an order of magnitude (3 of the eight 
GCs have inferred M/L < 0.1). This systematic under- 
estimate is not surprising, because the tidal forces expe- 
rienced by these halo globulars at this snapshot in time 
are smaller than the maximum tidal force that they have 
experienced in their past. Using their present distance 
in this M/L calculation will thus provide a lower limit 
on the masses required to yield the observed King tidal 
radii. Regardless, this calculation shows that Bootes II 
is an outlier from globular clusters with this metric of 
measurement. 

5. CONCLUSION 

In this paper, we use MMT/MegaCam imaging in g 
and r to present the first robust estimates of the fun- 
damental properties of the ultra-low luminosity Milky 
Way satellite Bootes II (d ~ 42 kpc). This object is old 
and its stellar population appears very similar to that 
of M92, showing that it is metal-poor ([Fe/H] < -2). 
With a total luminosity of only ~ 500 solar luminosities 

(M v 2.4±0.7 mag) and a half-light size of ~ 36±9 pc 

(assuming a Plummer profile), Bootes II lies away from 
globular clusters and dwarf spheroidals, but near Will- 
man 1 and Segue 1, in size-luminosity space. We showed 
that although the morphology of Bootes II appears ir- 
regular and elongated along the direction of the Galactic 
potential, that this distortion of its isodensity contours 
is not statistically significant in our dataset. 

The revised values we present for the distance, lumi- 
nosity, and physical size of Bootes II di f fer fro m those 
originally estimated in the IWalsh et al.l (|2007l ) discov- 
ery paper. The primary factor in these differences is 
the new distance estimate is 42 kpc rather than 60 kpc. 
The SDSS discovery data was more than two magnitudes 
shallower than the data presented in this paper. VLT 
photometry of Bootes II will be presented in Jerjen et al. 
(in prep.) along with a more detailed discussion of its 
stellar population and of the possible association, or lack 
thereof, with the Sagittarius Stream 



Our bootstrap analysis demonstrated the impact of 
small number statistics on the derived parameters for this 
ultra-faint class of objects, but showed that despite large 
uncertainties Bootes II has a size that makes it distinct 
from Milky Way GCs, although its King concentration 
is similar to that of Milky Way dwarf galaxies as well as 
to some diffuse GCs. 

The gap between Bootes II, Willman 1, and Segue 1 
and Milky Way globular clusters in size-luminosity space 
is not a selection effect because existing surveys would 
have been sensitive to such objects. However, the appar- 
ent separation in half-light size at 100 pc between dSphs 
and Bootes II, Willman 1, and Segue 1 could be a selec- 
tion effect, making it more likely that these three objects 
are fundamentally connected with the dwarf galaxy pop- 
ulation. We pointed out that old star clusters filling in 
that apparent gap may exist in environments other than 
the Milky Way, but would have escaped detection owing 
to their low surface brightnesses. 

If Bootes II is a self-bound system in equilibrium, it 
could represent the continuation of the dwarf galaxy 
population into the extreme low luminosity regime. We 
showed that it is reasonable to believe that Bootes II is a 
relaxed system, and used its King tidal radius to infer a 
lower limit M/L of ~ 98^4°. Dropping the assumption 
that the King tidal radius is physically meaningful for 
Bootes II, its spatial extent is larger than that naively 
expected for globular cluster-like mass-to-light ratios, if 
it is self-bound. However, we cannot rule out that it is a 
low M/L star cluster that is undergoing tidal disruption 
or disturbance, much like the Pal 5 cluster. 

Regardless of whether or not Bootes II is dark mat- 
ter dominated, it can provide a unique laboratory with 
which to investigate the low surface density limit of star 
formation and the tidal field of our Galaxy. The study 
in this paper also provides a partial road map for the 
future study of numerous similar objects that may be 
discovered in upcoming surveys f or Galactic s atellites to 
greater depth (e.g. PanSTARRS.lKaiser| [200l or in pre- 
vious ly unsearched sky (e.g. Skvmapper. iKeller et ail 
120071 and Walsh et al. in prep ). Ultimately, determining 
the boundedness of its stars will be needed to definitively 
pin down the nature of the peculiar Bootes II object. 
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TABLE 1 
Bootes II Properties 



Parameter Measured Uncertainty bootstrap median 



RA (h m s) 
Dec (d m s) 
(1,6) 

(to - M) 
Distance 
r h (PI) 
rh (ex) 
r h (PI) 
r h (ex) 
My (PI) 

M v (ex) 

Mo,v (PI) 
Mo,v (ex) 
r c 
n 
6 

e 
A 

A minor 



13 58 05.1 
+12 51 31 
(353.75,68.86) 

18.1 
42kpc 
2.8 arcmin 
2.5 arcmin 

35 pc 

31 pc 

-2.3 

-2.2 

27.76 mag arcsec - 2 
27.70 mag arcsec -2 

25 pc 

127 pc 

-33° 
0.27 



±7.2s 
±23// 

±0.06 
±1.6 kpc 
±0.7 
±0.8 
±9 
±10 
±0.7 
±0.7 
±0.31 
±0.33 
±9 
±35 
±57° 
±0.15 
±1.01 
±1.03 



13 58 04.3 
+12 51 09 

18.1 
42 kpc 
2.8 arcmin 
2.6 arcmin 

36 pc 

33 pc 

-2.4 

-2.2 

27.93 mag arcsec -2 
27.90 mag arcsec -2 

25 pc 

155 pc 

-28.5° 
0.34 

-0.26 

-0.33 
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